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ATMOSPHERIC  IR  TRANSMISSION  DATA  KOK  A  TEMPERATE  MARITIME  ENV I E'.AM.'.A  ; 

D-K.  Cutton 

S  C  M  M  A  H  Y 

"* Measurements  are  reported  on  infrared  tr.insfr.iYsMoj;  lErou^E 
the  atmosphere  at  a  (.oastal  marine  site  ii;  vi  i«'T:pt‘:Mti' 
environment.  Ten  broadband  spotitral  rc'^jons  wi-re  asr.j  j:-..: 
variations  of  tranjimittance  with  absolute  h'.imi.iiiy  «iud 
visibility  over  paths  of  5  and  km  invesi  M .  .Ml 

data  hav'e  been  compared  with  tht‘  pr**-!  u  i  lon'^  i;,.  AEMl. 

computer  code  EOWTR.AN  .  The  d  isi.reparu.y  ri'p.'* 
previously  between  tlje  measured  arel  caiuulaied  :  ii  i 
4.4  to  5.4  pm  region  still  exists,  .iltE.ough  :.-'u  <:  r.  »  «  •. 
the  water  continuum  absorption  coel  f  1 1,  i  rut  s  !->r 
region  have  reduced  this.  Reasonable  agreement  his  h,-. 

shown  to  exist  in  the  8  to  12  um  region  for  abs.  .i- 
humidities  up  to  12  gm/m’.  The  results  of  the  eompar  s 
for  the  1  to  2.5  pm  region  showed  a  greater  d:s^repi:<y. 

Some  of  this  discrepancy  coul<l  be  arising  :  ro:n  jK-.nt 
measurements  of  aerosol  extinction  using  imju  .pmctii  wl  ih 
is  underestimating  the  aerosol  dtt«'nuat  lou  r.oi- f !  i c  nun 
Data  on  scintillation  fluctuations  recorded  dining  the 
transmittance  measurements  are  also  presented. 
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I.  r.STKODUCTiO.V 

Data  on  broadband  infrared  (IR)  atmospheric  transmittance  have  been  colIicLi-d 
over  long  paths  near  the  air-sea  interface  in  a  temperate  environment  di.rmg 
the  period  November  1977  to  December  1979.  This  work  was  carried  out  under 
task  NAVY  77/141  which  involved  a  study  of  the  parameters  affecting  '.K 
atmospheric  transmi ttance  near  the  air-sea  interlace.  Such  data  in 
particular,  are  required  to  verify  the  predictions  of  atmospheric 
transmittance  models  e.g.  the  AFGL  computer  code  LOWTKAN  9,  under  similar 
meteorological  conditions  ( ref .  1 )  .  The  operating  c.ipab  1 1  ;  t  ics  of  Is 
surveillance,  missile  and  countermeasure  equipment  in  the  env  i  ror.monts 
encountered  around  Australia  are  roost  important  to  Australia's  defence  and  the 
evaluation  of  such  systems  generally  depend  on  the  use  of  atmospl'.e  r  ic 
transmittance  models  to  predict  the  optical  properties  of  the  intervening 
atmosphere.  Hence  the  need  exists  to  validate  the  predictions  of  models  over 
long  ranges  in  Australian  environments. 

A  preliminary  assessment  of  transmittance  data  me.isiired  ,u  the  coastal  .mai  i.-.e 
site  over  the  period  November  1977  to  April  1978  lius  already  hee:. 
reportedfref .  2) .  It  is  the  purpose  of  this  paper  to  report  on  the  an  wy.sis  of 
transmittance  data  measured  during  the  period  October  1976  to  Decemuer  l')7‘j. 
Data  are  included  for  ten  broadband  spectral  regions  from  I  to  12  pm  anu  the 
effects  of  absolute  humidity,  relative  humtlity  and  visibility  on  l.s 
atmospheric  transmittance  are  examined.  Predictions  from  LOWT.RAN  i  tor 
similar  meteorological  conditions  are  made  and  the  corresponding  d.it.; 
compared.  A  brief  qualitative  analysis  is  .also  reported  on  some  scint  i  i  1  it  ;  o:i 
data  recorded  during  IR  transmission  measurements. 


2.  INSTRLMENTATION 

The  instrumentation  used  to  measure  absolute,  broadband  IK  atmcsphi  :- . 
transmittance  and  basic  meteorological  par.arooters  has  already  beer,  descr ;  is  c. 
in  detail  (ref .  2 , 3) .  In  brief,  the  transmission  measuring  equipment 
of  a  trailer-mounted  broadband  IR  source  and  a  van-mounted  radiometer  ab'..-  •,  o 
operate  up  to  ranges  of  10  km.  The  system  can  mcasuii'  the  absolu;  e 
atmospheric  transmittance  in  selected  spectral  regions  using  a  combination  of 
three  cooled  detectors  (PbS,  InSb  and  MCT)  and  spectral  filters.  The 
efficient  IR  source  is  a  re-entrant  conical  black-body  cavity  with  an  jpe.-iun 
of  100  cm’,  a  nominal  surface  temperature  of  975  K  and  a  lota!  radiant  power 
output  of  approximately  500  W.  Both  the  !R  source  and  radiometer  are 
calibrated  as  described  in  reference  3.  Appendix  I  gives  the  derivation  of 
the  transmissometer  measurement  equation  while  a  summary  of  the  error  analysis 
performed  on  this  system  is  given  in  Appendix  II. 

During  the  transmittance  measurements  the  air  temperature,  relative  humidity, 
visibility,  sea  state,  wind  speed  and  wind  direction  were  monitored.  Wet  and 
dry  temperatures,  and  aerosol  extinction  coefficient  (at  X  =  0.55  pro',  were 
measured  at  the  two  ends  of  the  transmission  path  at  10  min  intervals  and  the 
data  averaged. 

Atmospheric  particulate  matter  was  sampled  during  measurements  at  the 
radiometer  end  of  the  path  using 

(i)  47  mm  diameter  Miliipore  filters  with  0.1  pm  pore  size  at  flow  rates 
of  approximately  4.5  L/min  and 

(ii)  120  mm  diameter  Whatman  41  filter.s  with  flow  rato.s  of  approximately 
0.5  m’ /min . 
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kosults  from  lIk'  analysis  of  tlioso  filters  uiil  In*  liu*  sulijiu.i  o;  l.iii-i 
reports.  Briefly,  the  samples  on  the  Milliporo  filters  provuiod  aerosol  s-./.r 
distribution  curves  and  the  Whatman  41  samples  uore  am.  lysed  usin^  \-r.iy 
fluorescence  techniques  to  determine  the  elemental  composition  of  the  aerosol, 
quantitatively. 


3.  SITE  DESCRIPTION 

The  atmospheric  transmittance  and  meteorological  meaj>urcments  were  made  at  t:io 
same  coastal  site  referred  to  in  reference  2.  This  site  was  chosen  as  best 
suiting  the  requirements  for  open,  occan-type  sea  conditions,  ie  it  has 
prevailing  on-shore  winds  not  passing  immediately  over  continental  land  masses 
and  has  a  low  industrial  contamination.  Figure  1  shows  tlie  relative  location 
of  the  site  with  respect  to  Adelaide,  while  figure  2  gives  details  of  the  two 
transmission  paths  used  (5.03  and  9.05km)  and  the  topography  of  the  site. 


4.  DATA  ANALYSIS 

The  data  reported  in  this  paper  cover  the  period  from  October  1978  lo 
December  1979.  Transmittance  measurements  were  made  on  o2  days  during  this 
period  for  5.03  and  9.05  km  pathlongths  and  in  the  spectral  regions  defined  by 
the  filter  half  power  points  listed  in  Table  1  below. 


TABLE  1.  BROADBAND  SPECTRAL  REGIONS  SELECTED  FOR  ATMOSPHERIC  TRANS.MITT.\NCE 
MEASUREMENTS 


Filter  No. 

Half  Power 

Points  (pm) 

F i Iter  No. 

- 1 

Halt  Power 

Po i nts  ( wm  ) 

1 

1.A80-2..300 

8 

4.410-  .3.410 

2 

1 .4A0-1 .820 

10 

8. 200-11. 7n0 

3 

1 .910-2.470 

11 

8.330-  9.800 

5 

3.5S0-4.000 

12 

10. ..70-10.  700 

7 

4 . 3b0-4 . 590 

13 

0.9.33-  1.150 

Figure  3  reproduces  the  spectral  curve.s  of  the  filters  useii. 

All  transmission  and  meteorological  data  were  stored  on  a  IBM  370/3033 
computer  and  programmes  were  written  to  sort  and  plot  the  transmission  data 
into  the  required  parameter  groupings.  The  transmittance  data  were  derived 
from  the  recorded  radiometer  output  signals  which  spanned  at  least  a  3  min 
interval  and  which  were  averaged  by  eye.  Measurements  for  each, 
filter/detector  combination  were  not  made  at  regular  intervals  but  generally 
repeated  as  soon  as  the  filter  wheel  completed  one  revolution.  (The  ma.ximum 
number  of  filters  which  could  be  held  in  the  wheel  at  one  time  was  seven). 
The  radiometer  data  were  converted  using  the  appropriate  calibration  factors 
measured  in  the  laboratory  for  each  f i 1 tcr/detector  combination  used.  Those 
data  where  there  was  an  indication  that  precipitation  possibly  existed  in  the 
path  were  not  included  in  the  data  set. 

The  data  were  sorted  to  show  the  variability  of  atmospheric  transmittance  with 
(i)  water  vapour  content  for  filters  5  to  12  and 


(ii)  visibility  for  filters  1,2, 3, 5  and  13, 

Further  subdivision  of  the  data  group  related  to  water  vapour  content  was  not 
undertaken  because  of  the  limited  range  of  air  temperatures  measured  (13  to 


« 


f 


f 


3 


ERL-026  5 -T« 


/' 


23®C)  and  the  insufficient  number  of  data  when  smaller  temportituro  inturvais 
Were  chosen  (see  Section  S.l).  The  data  related  to  visibility  was  furtlu'r 
subdivided  to  accept  only  those  readings  associated  with  on-shore  winds  which 
at  this  Site  generally  fell  within  SSW  and  ESE.  Because  only  limited 
variations  in  wind  speed  (0.5  to  10  ms)  and  sea  state  conditions  wore 
experienced  no  division  of  the  data  was  made  with  these  parameters. 

It  is  interesting  to  note  that  visibilities  (=  3.91/0  ,)  as  measured  with. 

aerosol 

integrating  nephe lometers  at  two  locations  along  the  siglit  line  wore  greater 
than  50  km  on  those  days  measurements  were  made.  This  is  not  unexpected  as 
the  site  is  located  where  fogs  and  low  maritime  visibilities  are  not  generally 
experienced.  The  temperature  and  absolute  humidity  were  monitored  for  each 
end  of  the  path  and  averaged  as  described  in  reference  2.  The  average  wir-d 
angle  over  a  particular  hour  was  allocated  one  of  eight  45®  intervals  while 
the  corresponding  average  wind  speed  was  assigned  a  Beaufort  number.  Tliese 
classifications  are  summarised  in  Table  2. 


TABLE  2.  CLASSinCATlON  OF  WIND  SPEED  AND  DIRECT  ION 


Wind 

Classification  No. 

Wind  Angle 
Interval 
(degrees) 

- 1 

Wind  S-piM-d  1 

Interval  j 

i  m.s  j  ! 

1 

338-  23 

0 . 3  -  1  .  3  ' 

2 

23-  o8 

1  .  b  •  3  .  '3 

3 

68-113 

3  .  -» •  5  ,  *• 

U 

113-138 

3  .  3-  ■  .  '> 

5 

156-303 

8 . - 10 , : 

6 

203-248 

1  ,  b- : 

7 

248-203 

13,0-;-.'. 

8 

-  —  j 

293-338 

Figures  4  and  5  reproduce  plots  of  transmittance  versus  ahsoluie  l.uiti  i.iy 
content  covering  all  visibilities  encountered  for  six  .spectral  regions  n. iwc 
pathlengths  studied.  Transmittance  data  shown  in  the  two  i  iguros  w^rc, 

however,  grouped  into  four  temperature  intervals.  Figures  n  to  8  .•"eproueu 
plots  of  transmittance  versus  visibility  where  the  data  for  each  group 
been  divided  into  five  relative  humidity  intervals. 

In  examining  more  closely  the  effects  of  water  vapour  content  on  transmitta:...* 
a  similar  procedure  was  used  as  described  in  reference  2.  Transmitt.ince  date 
for  a  particular  spectral  region  and  visibility  grouping  were  divided  in’o 
equal  intervals  of  0.5  gm/m’  absolute  humidity  and  the  data  falling  wii!..;. 
each  interval  averaged.  Hence,  in  figures  9  to  11  and  13  to  15  the  v^rossos 
represent  the  average  values  and  the  bars  the  extreme  data  limits  whilst  ihe 
number  of  data  falling  in  each  interval  are  given  in  the  figure.  The  absolute 
humidity  interval  of  0.5  gm/m*  was  selected  as  it  represents  the  maximum 

uncertainty  likely  to  be  experienced  in  the  measurement  of  absolute  humidity 
along  the  path. 

For  the  detailed  analysis  of  the  effects  of  the  maritime  aerosol  on 

transmittance  in  the  1.0  to  4,0  ym  spectral  region,  the  data  for  each  range 
were  grouped  into  equal  visibility  intervals  of  20  km  and  five  equal  intervals 
of  relative  humidity  between  45  and  95%.  Data  falling  in  each  visibility  and 
relative  humidity  interval  were  averaged  and  figures  19  to  2b  reproduce  these 
averaged  data.  The  crosses  represent  the  average,  values  and  the  circles  the 
extreme  data  limits.  The  contribution  of  the  aerosol  extinction  coefficient 
to  the  total  extinction  coefficient  in  the  4.4  to  S.4  yro  and  8  to  12  urn 
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spi-ctral  rogion  is  small  compared  to  the  molecular  extinction  for  those 
visibilities  experienced  at  the  V'ictor  Harbor  site.  The  lower  visibil  ty  .i: 
SO  km  recorded  at  the  site  gives  maritime  aerosol  extinction  coefficients  o'. 

approximately  0.020  km  '  and  0.014  km  *  for  the  two  spectral  regions  at  .i  SOT 
relative  humidity,  as  compared  with  the  measured  lower  limit  of  effective 

total  extinction  coefficients  of  0.3  and  0.1  km  respectively.  Tor  '.hese 
reasons  no  plots  of  transmittance  versus  visibility  are  provided  for  these  twe 
regions  and  as  well  no  distinction  was  made  between  off-shore  and  or;-sl-.ore 
winds  for  the  visibility  groups  in  figures  0  to  15. 

4.1  LOVi’TRAN  .5  computations 

The  analysis  has  included  the  comparison  of  measured  transmittance  with  li.e 
predicted  values  from  LOWTRAN  5  for  similar  meteorological  conditions.  In 
this  analysis  the  data  have  been  presented  in  both  the  form  o;  effective 
transmittance  and  effective  extinction  coefficients  for  a  parti  .'.ii 
wavelength  interval.  The  effective  traiismi  ttance  computations  using 
LOWTRAN  5  were  carried  out  by  the  numerical  integration  of  the  iv 
integrals  given  in  equation  (1.5)  in  Appendix  I;  these  dat.i  are  represi'iiteu 
by  the  solid  lines  labelled  LT5  in  the  figures.  Exceptions  to  tins  geiiei-.i'. 
rule  are  figures  14  and  15,  whore  LT5A  computations  are  represented  by  ,> 
solid  line  and  LT5  computations  by  a  dashed  line.  Plots  of  measureu  ter-  ;s 
calculated  effective  extinction  coefficients  are  also  presented.  These 
plots  provide  a  possible  means  of  determining  the  deviation  of  the  s'.oiie  of 
a  least -squares  fitted  curve  from  the  lil  situation;  the  fitted  curves  .lie 
represented  by  the  dashed  lines  in  figures  12  and  16.  The  cor.tr  ibu;  ;c:. 
from  the  regression  of  the  calculated  extinction  coefficients  cn  t:.> 
measured  coefficients  (coefficient  of  delorminat  ion  r’ )  can  also  be 
computed.  It  should  be  pointed  out  that  linear  regression  analysis  is  no', 
strictly  applicable  when  uncertainty  in  both  sets  of  data  is  compar.ihle  bn; 
if  caution  is  used  in  interpreting  the  results  some  indication  of  trends 
can  be  ascertained.  (The  EOWTRAS  code  has  been  di-rived  to  model  low- 
resolution  atmospheric  transmittance  to  within  lOT) .  It  was  furtler 
assumed  that  the  fitted  linear  curve  passed  through  the  origin.  This  would 
be  true  if  the  calibration  constants  for  each  f i Iter/dctector  corabinatio: 
were  correct,  and  the  trend  of  the  data  in  the  plots  indicated  that  tins 
was  generally  close  to  being  the  case,  particularly  in  the  8  to  12  tni 
region  where  smaller  extinction  coefficients  arc  experienced. 

It  should  be  noted  that  the  effective  extinction  coefficient  is  the  sum  oi 
the  aerosol  and  effective  molecular  extinction  coefficients.  The  eftective 
molecular  extinction  coefficient  is  determined  for  a  particular  patliieiigth 
and  wavelength  interval  AX.  As  many  direct  spectral  absorption  lir.es 
occurring  within  AX  are  not  resolved  the  Bougor  exponential  transraittar.ee 
law  =  exp(-oR)  cannot  generally  be  used.  This  arises  because  the  fine 

structure  of  the  spectrum  varies  differently  with  pathlcngth  as  a  result  of 
selective  absorption,  and  thus  the  effective  molecular  extinciior 
coefficient  cannot  be  used  to  compute  the  molecular  atmospheric 
transmittance  for  other  pathlenglhs  and  wavelength  intervals.  The  aerosol 
extinction  coefficient,  however,  only  varies  slowly  with  wavelength  and 
docs  obey  the  exponential  transmittance  law. 


5.  RESULTS  AND  DISCUSSION 

This  section  will  describe  the  results  of  the  analysis  of  variation  ot 
transmittance  with  absolute  humidity  and  visibility.  One  general  comment  can 
be  made  about  the  scatter  that  is  always  evident  when  atmospheric 
transmittance  measurements  are  made.  It  is  readily  apparent  that  the  scatter 
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in  the  reduced  transmittance  data  for  the  October  1978  to  December  1979  period 
is  less  than  the  scatter  in  the  corresponding  data  measured  between 
November  1977  and  April  1978  and  reported  in  reference  c.  During  a  laboratory 
calibration  carried  out  between  these  two  periods  it  was  found  that  ti.i' 
alignment  of  the  detectors  in  the  radiometer  was  not  correct  and  this  appeared 
to  have  resulted  in  signals  from  the  source  not  always  being  correctly  peased. 

5.1  Variation  of  atmospheric  transmittance  in  the  8  to  12  urn  region 

In  this  section  the  results  obtained  with  filters  10.11  and  12  are 
discussed.  The  meteorological  parameter  which  has  most  influence  on 
atmospheric  transmittance  in  this  spectral  region  is  the  absolute  humidity 
or  the  water  vapour  content  in  the  path.  Although  this  region  does  not 
contain  many  strong  spectral  absorption  lines  arising  from,  the  uniform  gas 
mixture  and  water  vapour,  there  is  a  strong  anomalous  absorption  whi,,h 
varies  only  slowly  with  wavelength,  and  which  dominates  the  transmitt.ir.ee. 

Modelling  of  this  anomalous  absorption  by  water  vapour  in  LOWTKAS  5  is 
based  on  an  expression  for  water  continuum  absorption  derived  empir  i,;.i '.  ly 
from  laboratory  measurements  on  water  vapour  absorption  in  short  path '.  e.ugth 
cells.  This  expression,  summarised  in  reference  A,  exhibits  a  qujir.itic 
dependence  on  water  vapour  pressure.  However,  the  fundamental  form  tr.  is 
molecular  absorption  mechanism  takes  is  still  subject  to  conjecture,  and 
currently  two  theories  exist  about  the  mechanisms  invo I vedfre f . 5 , n  , 

Laboratory  and  field  measurements  to  date  have  indicated  that  the  omp  i  r  ic.t '. 
model  used  in  LOWTRAN  5  predicts  the  extinction  coefficient  for  anoma '.ous 
water  vapour  absorption  to  a  reasonable  degree  of  accuracy  for  w  iter  v,i:..,  .r 
content  up  to  15  gm/m’  for  air  temperatures  in  the  region  of  2C  to  .. 

There  is  not  sufficient  data  with  a  wide  spread  of  temperature  to  .itie...-;  , 
with  accuracy,  the  negative  temperature  (exponential  Solta.m.in:,  - le  ■ 
dependence  of  the  form  exp{A(296/T-l))  that  this  anom.a'.oiis  ahsorpi  .... 
been  predicted  to  follow. 

The  transmittance  data  obtained  from  the  Victor  Harbor  site  are  reprodu,.eu 
in  figures  4(a),  4(b)  and  4(c)  as  plots  against  absolute  humidity, 
these  plots,  the  transmittance  data  were  further  subdivided  .;.ic 
temperature  intervals  of  5’C.  It  is  evident  that  most  of  the  data  are 
associated  with  air  temperatures  ranging  from  13  to  23°C.  In  view  of  t.te 
small  temperature  range  and  the  small  number  of  data  falling  into  a  given 
temperature  and  absolute  humidity  interval,  it  was  decided  not  to  perform 
any  analysis  with  temperature  as  a  parameter.  Computations  undertaken  wii;. 
the  current  version  of  LOWTRAN  provide  transmittance  values  that  change  ny 
a  maximum  value  of  0.04  transmittance  at  the  maximum  measured  absolute 
humidity  of  15  gm/m’  in  the  temperature  range  13  to  23°C.  Hence,  it  woula 
be  unlikely  such  changes  in  transmittance  with  temperature  could  be 
detected  in  the  present  data  unless  a  large  number  of  measurements  were 
taken  to  reduce  to  a  minimum,  the  scatter  in  transmittance. 

The  plots  of  averaged  transmittance  data  for  each  equal  absolute  humidity 
interval  (figures  9  to  11)  indicate  reasonable  agreement  with  LOWTRAN  i 
predictions.  To  quantify  these  differences,  plots  of  calculated  versus 

measured  effective  extinction  coefficient  (refer  to  Section  4.1)  were  done. 

Figures  12(a),  12(b)  and  12(c)  reproduce  those  plots  whore  the  a.irosol 

extinction  coefficient  was  less  than  0.04  km  This  figure  was  chosen 

because  most  of  the  visibilities  occurred  in  this  region  and  since  the 
aerosol  extinction  coefficient  is  small  (even  for  relative  humidities 
around  80%)  it  would  only  contribute  a  small  uncertainty  to  the  I 

transmittance.  The  results  of  the  linear  regression  analysis  are 
I  summarised  in  Table  3.  , 

The  values  if  b,  the  lope  of  the  regression  line,  indicate  that  the 
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TABLE  3.  SUKS^ARY  OF  LINEAR  REGRESSION  ANALYSIS  FOR  THE  8  TO  12  Mm  SPECTRAL 
REGION 


(a)  Visibility  <  96  km 


Filter 

No. 

Spectral 

Region 

(Mm) 

Range 

5.03  km 

9.05  km 

b 

r^ 

SD 

Uncertainty 

in  T  „ 
eff 

- 1 

b 

r^ 

SD 

Uncertainty 

in  T  j.. 
eff 

10 

8.20-11.76 

1.12 

0.997 

0.009 

+0.05 

-0.09 

1.12 

0.975 

0.028 

+  0.20 
-0.  12 

11 

8.33-  9.80 

0.977 

0.995 

0.012 

±0.06 

1.03 

- 1 

0.9"5 

0.026 

— 

+  0.  IS 
-0.11 

12 

10.47-10.70 

0.921 

_ 

0.993 

_ 

0.015 

+0.08 

-0.07 

1.07 

_ 

0.974 

_ 

0.029 

+  0.21 
-0.12 

(b)  Visibility  >  96  km 


Filter 

No. 

Spectral 

Region 

(Mm) 

Range 

5.03  km 

9.05  km 

b 

r' 

SD 

Uncertainty 
in  T 

eff 

b 

r' 

SD 

Uncertainty 

in  T  . 
ett 

10 

8.20-11.76 

1.20 

0.997 

0.01 

+0.06 

-0.04 

1. 13 

0.992 

0.015 

+0.09  1 

-0.07 

11 

8.33-  9.80 

1.12 

0.982 

0.02 

+0. 12 
-0.08 

1.05 

0.995 

0.011 

1 

±0.00  1 

12 

10.47-10.70 

0.999 

0.985 

0.02 

+0.12 

-0.08 

1.09 

0.983 

0.02 

+0.1.  . 

-0.09  1 

average  fit  of  the  calculated  data  to  the  measured  data  is  within  ' 
except  for  one  filter/range  combination.  Assuming  an  uncertainty  of  plus 
or  minus  two  standard  deviations  between  the  calculated  and  measured 

extinctions  and  assuming  mean  transmittances  of  0.5  and  0.3  for  the  two 

ranges  5.03  and  9.05  km  respectively,  then  the  uncertainty  in  the  effective 
extinction  coefficient  can  be  translated  to  an  uncertainty  in  the 

transmittance.  The  values  listed  in  Table  3  reveal  a  reasonably  consistent 
uncertainty  in  transmittance  for  the  5.03  km  range  while  there  appears  to 
be  a  higher  uncertainty  for  the  9.05  km  range. 

5.2  Variation  of  atmospheric  transmittance  in  the  3.5  to  5.5  pm  region 

The  results  discussed  in  this  section  cover  the  data  obtained  for 

filters  5,7  and  8  (see  figures  5(a),  5(b)  and  5(c)).  The  main  meteor¬ 
ological  parameter  influencing  the  transmittance  is  again  water  vapour  with 
maritime  aerosol  extinction  also  contributing  but  to  a  lesser  extent.  This 
region  is  characterised  by  strong  absorption  due  to  the  many  spectral 
absorption  lines  arising  from  HiO  and  CO,  molecules  as  well  as  a  N, 

continuum  absorption,  and  an  anomalous  absorption  due  to  water  vapour 
continuum. 


5.2.1  -*.4  Lo  5.4  Jim  spnctral  region 


Referring  to  pret'ious  investigatioiib  by  this  author  in  the  to  5 .  -  -t. 

spectral  region(ref  .4)  it  was  found  that  the  addition  of  water  v.i:v'..r 
continuum  absorption  to  the  LOKTRAN  model  improved  con.s  iderab  1  y 
accuracy  of  the  predictions.  Bcn-Shalom  et  al  in  a  recent  paper , vef . 
also  reported  spectral  transmittance  measurements  with  a  resolui'.cr.  oi 
0.07  ym  for  the  4.3  to  5.5  ym  region  which  showed  a  similar  disc:'i  p..i:. -y 
when  compared  with  computed  transmittances  from  LOWTKAS  4.  .M'tor 
included  the  water  continuum  absorption  o.xprcss  ion  into  1>0mTS..\N .  j.s 
suggested  by  this  author(rcf  .4) ,  good  agreement  is  obtair.ed  for 

pathlenglhs  near  b  km  although  it  is  not  clear  whether  they  extended  liie 
water  continuum  absorption  expression  down  to  4.3  ym.  Tin*::- 
experimental  measurcmentsiref  .fe)  also  provided  values  of  t;;e  srlt- 
broadening  absorption  coefficient  in  the  water  continuum  aosorptio.. 
expression  which  were  in  agreement  (to  within  about  SO^o)  with  those 
derived  by  this  author.  Considering  the  various  experimental  errors 
involved,  this  agreement  is  encouraging. 

However  the  transmittance  data  reported  here  for  tlic  4.**  to  5.-  -.r. 
region  (figures  14  and  15)  indicate  that  the  measured  effect. v»' 
transmittance  is  lower  than  that  given  by  ca’culation  using  hOUT^.VN  t 
(see  LT5  curves  in  the  figures)  and  the  discrepancies  are  depe..d.  :.t 
pathlength  and  absolute  humidity.  It  would  appear  th.at  the 
continuum  absorption  model  included  in  LOWTRA.S  5  may  not  .i.-..ou:.t 
completely  for  the  discrepancies  observed  in  this  spectral  region. 

Recent  ly  ( ref  .9)  AFGL  have  examined  this  region  and  »..omputi:d  .. 

self  and  foreign  broadening  absorption  coof  f  Ic  ients  i,C^(v,7- 

Cp.(v,T))  over  the  region  from  v=i00  to  3000  cm  *  duo  to  :ur 
absorption  of  HjO  absorption  lines.  Tiieso  curves  have  b»?eii  wO 

compute  the  water  continuum  transmittance  curve  in  figure  17(a')  for  .';.o 

conditions  stated,  over  1800  to  3000  cm  ^  region.  Also  included  iT' 
transmittance  curves  derived  from  the  water  continiium  abscrci  ur. 

coefficients  used  in  *^0WTRAN  S  for  2400  to  3000  cm  ^  region,  and  ny  ihi' 

author  for  the  2000  to  2200  cm  '  region.  In  view  that  the  data  in  thf 

2000  to  2200  cm  '  region  were  derived  from  field  measurements  the 
agreement  is  reasonable.  Figure  17(b)  shows  a  plot  of  Y  = 

V  derived  from  the  calculated  C^(v,T)  and  Cj_.(v,T)  in  reference  9.  It  ;s 

noted  that  Y  varies  over  the  range  of  0.02  to  0.12  and  it  is  clearly 
evident  that  water  continuum  calculations  for  this  region  should  inclun* 
F  as  a  function  of  v.  (LOWTKAN  5  assumes  a  fixed  value  of  0.12  for  tin* 

2400  to  3000  cm  '  region.)  In  figures  14  and  15  the  LT5  curves  refer  to 

computed  curves  using  LOVTRAN  5,  and  LT5A  curves  to  LOkTRAN  5 
incorporating  the  AFGL  water  continuum  absorption  data  over 

1800  to  3000  cm  ^  region  and  J  varying  with  wavelength  as  in 

figure  17(b).  The  maritime  aerosol  model  was  used  for  all  LOWTRAN 
coroputat ions . 

In  the  case  of  filter  8  the  agreement  of  the  data  arc  reasonable  while 
for  filter  7  there  appears  to  be  evidence  of  a  significant  discrepancy 
remaining.  In  order  to  quantify  the  differences  between  the  measured 
and  calculated  data  using  the  new  continuum  data,  a  linear  regression 
analysis  was  performed  on  the  calculated  and  measured  effective 
extinction  coefficients  for  filters  5,  7  and  8.  (Plots  of  these  data 
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TABLE  4.  SUMMARY  OF  LINEAR  REGRESSION  ANALYSIS  FOR  THE  3.5  TO  5.5  vtr,  REOION 


(a)  Visibility  <96  km 


F  liter 
No. 

Spectral 

Region(iim) 

Range  j 

5.03  km 

0,05  km  i 

b 

r’ 

SD 

b 

- 1 

- 

SD 

5 

3.55-4.00 

0.74 

0.978 

0.015 

0.7’ 

0.9')2 

0  .  ' 

7 

4.36-4.59 

0.84 

00.998 

0.018 

0.83 

0  .  'J')") 

0.00  5 

8 

4.41-5.41 

0.9t>5 

0.996 

0.018 

0.93 

0  .  ‘)'rj 

; 

_ 1 

(b)  Visibility  >9t>  km 


Filter 
No . 

Spectral 

Region(ijm) 

Range 

5.03  km 

9 . 0  5  km  i 

b 

SD 

b 

1 

r'  SO  ; 

5 

3 . 55-4,00 

0.95 

0.950 

0.018 

0.86 

0.99o  I  (.  .00  5  j 

7 

4,36-4,59 

0  .89 

0.996 

0.019 

0.81 

D.'ViS  j  1 

8 

4.41-5.41 

1.04 

0.99  7 

_ 1 

0 . 020 

0.93 

_ 

^'.9  T)  1  . 

_ -.1  _ i 

NB:  Effective  extinction  coefficients  calculated  u  ;  l';:  :.ew  Al'iL.  u  : 
continuum  absorption  dataCref.O) 

for  the  high  visibility  groups  are  shown  in  figures  loib)  ai:d  Ir  ...  . 
In  all  except  one  case  the  slope  b  of  the  fitted  curves  .ill  W'  re  '.ess 
than  1.0  (see  Table  4), 

Because  filter  7  has  a  narrowr  bandpass,  errors  m.iy  arise  tr..r.  t  :.e 
filter's  spectral  curve  which  could  introduce  a  sm.ill  constant  error 
into  the  calibration  constants  and.  combined  with  low  transmittance 
values,  cause  the  discrepancy  observed  in  figure  14.  lienee,  t;.e 
differences  in  computed  and  measured  t  ransmi  t  tance  for  the  s...  to  4 ,  n  urn 
region  were  examined  more  closely.  The  transmittance  dat.i  in  figure  I- 
were  fitted  with  linear  curves  using  linear  regression  analysis  and  tlie 
difference  in  the  calculated  and  measured  effective  extinction 
coefficients  determined  for  a  given  absolute  humidity.  The  results  ire 
given  in  Table  5.  Three  out  of  the  four  sets  of  data  exhibit  a  near 
constant  difference  with  water  vapour  in  support  of  the  argument  tuat  a 
constant  error  may  explain  the  discrepancy.  (It  is  pointed  out  that  the 
computed  transmittance  values  using  filter  7  reported  in  n  fereiice  2  are 
in  error  since  transmittance  was  calculated  between  wavelengths 
corresponding  to  the  filter  half  power  points  and  not  for  wavelengt!;s 
covering  the  whole  filter  spectral  range  as  now  done).  Any  discrepancy 
arising  in  the  LOWTRAN  computation  of  the  spectral  absorption  due  to  the 
band  model  used  for  the  uniformly  mi.xcd  gases  has  been  ruled  out 
following  checks  in  this  region  with  the  1978  AFGI.  atmospheric  line 
compilation  tape  HITRAN(ref .2) . 
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TABLK  5.  ESTIMATED  DISCREPANCY  IN  CAIXLL<VrED  AND  MEASURED  EFFECTIVE 

EXTINCTION  COEFFICIENTS  AS  A  FUNCTION  OF  ABSOLUTE  HUMIDITY  FOR 
4,^  TO  4.6  pm  REGION 

(a)  Visibility  <96  km 


Range 

Absolut 

0  humidity  | 

(km) 

6 

9 

12 

5.03 

0.072  km'* 

0.077  km  ^ 

0.077  km'* 

9.05 

0.053 

0 . 

0.061 

(h)  Visibility  '9n  km 


Range 

Absolute  humidity  i 

gm  m*  )  1 

(km) 

b 

9 

i  2 

5.03 

0.03b  km  ^ 

0.048  km  ^ 

Q.O 58  km  '  ; 

9.05 

0 . 030 

O.Ool 

0 .  i 34  1 

,  J 

5.2.2  3.5  to  4.2  pm  spoctr.il  rollon 

For  the  transmittance  d.it.a  moasiirod  in  tbo  3.5s  to  - . u.m  ro*-..>:’. 
regression  analysis  reveals  a  greater  v.iriability  in  the  slope 
fitted  curves  (see  Table  4  and  figure  lo(a)).  Some  of  this  vari.it-.c.. 
must  arise  from  experimental  errors  when  t r.insm i ttan., e  values  are 
(see  Appendix  II).  For  water  vapour  concentrations  up  ’.  .li.-.  .t 

12  gm/m*  ,  the  Haught  and  Cordray  and  Bcn-Shalom  et  al  spect.  U. 
transmittance  data ( ref  s  .  7  , 10  and  11)  measured  over  5  to  o  km  path  lengt  ..s 
have  shown  reasonable  agreement  with  LOWTRAN  5.  There  is  evidence, 
however,  of  a  discrepancy  in  both  groups  of  data  when  water  vapour 
concentrations  exceed  12  gm/ra’ .  In  reference  4  it  was  reported  tiiat 
there  existed  a  discrepancy  between  laboratory  measurements  of  ti.e 
3.5  to  4.0pm  H2O  continuum  absorption  coefficients  at  Uvo  K  ny 

White  et  al(ref.l2)  and  Burch  ct  al(ref.l3). 

More  recently,  Watkins  ct  al(ref.l4)  carried  out  further  laboratory 
measurements  (using  the  Atmospheric  Sciences  Laboratory  21  m  absorpticr. 
cell)  to  determine  more  precisely  the  nature  of  the  continuum  absorplie:. 
at  ambient  temperatures.  Their  results  indicated  that  there  are 
contributions  from  both  far  wing  and  aggregate  water  molecule 
(eg  water  dimer  or  ion  cluster)  types  of  absorption.  In  contrast  the 
Burch  data  for  296  K  (on  which  the  continuum  absorption  model  used  in 
LOWTRAN  is  based)  arc  derived  by  extrapolation  from  data  measured  at 
338  K  where  apparently  the  predominate  absorption  is  due  only  to  the  tar 
wing  type,  ie  the  LOWTRAN  model  under-estimates  the  tranbmittai],.e 
because  of  the  missing  aggregate  water  molecule  type  of  absorption.  The 
inclusion  of  this  type  of  absorption  makes  the  overall  continuum 
absorption  more  strongly  dependent  on  water  vapour  pressure  than  is 
presently  shown  by  the  model  in  LOWTRAN.  As  a  result  of  the  reduced 
dependence  of  the  water  vapour  continuum  absorption  on  foreign  pressure 
broadening  the  absorption  will  not  fall  off  as  rapidly  with  combined 
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Nj-O,  pressure  (eg  with  increasing  altitude  as  is  the  case  with  tin- 

present  model).  As  the  two  types  of  absorption  contribute  appro.x  im  ate y 
equal  proportions  at  24°C  the  continuum  extinction  coefficients  are  now 
double  those  given  by  the  model  in  LOUTRAN'. 

To  sec  whether  this  new  model  will  go  some  way  to  .iccocnt ; t.g  for  t:.o 
discrepancy  reported  in  the  overseas  data  measured  at  iiigi.  w.it.-r  vapo..;- 
pressures  an  illustrative  example  is  provided.  Thi-  dat.i  from  T.ii  iv  o  c; 
Watkins  et  al  for  80“o  relative  humidity  (whicli  corresponds  to 
17.25  gm/m’  at  24°C)  were  used  to  compute  a  correction  factor  it 
selected  wavelengths  to  the  continuum  absorption  data  used  in  LuWTRAS. 
These  factors  were  applied  to  the  5.1  km  pathlengtii  LOWTKAS  curve  give;, 
in  figure  1(c)  of  Haught  and  Cordray.  The  modified  LOWTKAS  data  points 
are  reproduced  in  figure  18  along  with  the  spectral  transmittance  curve 
originally  measured  by  Haught  and  Cordray.  As  no  mention  is  made  of  the 
temperature  at  which  the  original  LOWTKAS  computations  wore  made,  the 

temperature  of  2a‘’C  had  to  bo  used.  Except  (or  the  2700  to  2800  cm  ‘ 
(3.70  to  3.57  pm)  region  the  transmittance  predictions  of  the  modifie.l 
LOWTRAS  model  are  closer  to  the  measured  values. 

The  discrepancy  in  the  2700  to  2800  cm  *  region  is  not  supported  by  the 
HjO  continuum  transmittance  data  in  figure  17ia)  based  on  the  new  AKGL 

computed  continuum  absorption  data.  Data  given  in  that  figure  for 

2600  to  2700  cm  '  region  does  indicate  about  a  O.Oj  correction  iacto: 
needs  to  be  applied  to  the  present  LOWTRAS  5  HjO  data. 

Furthermore,  the  temperature  dependence  of  the  HjO  continuum  ahsorp.t  ; 

for  this  spectral  region  requires  more  laboratory  i nvest  ig.it  ions  : 
define  the  precise  relationship.  The  Watkins  et  al  HjO  continuum  model 

predicts  that  the  absorption  will  sliow  a  much  stronger  d.epen.ience  cn 
temperature  near  20°C  than  the  Burch  data  now  presently  used  in  L.IWTKAS . 
Also  insufficient  field  data  on  transmittance  exists  in  this  regio:.  to 
examine  ttie  temperature  dependence  of  HjO  absorption.  The  present  d.it.i 

are  limited  in  number  and  too  largely  scattered  for  doing  any  .in.ily^is 
of  transmittance  variability  with  temperature. 

It  is  interesting  to  note  that  Haught  and  Cordray  derived  an  a  1 1  e  r::.r.  ;  ve 
algorithm  by  fitting  a  polynomial  to  their  data  for  this  spec'.r.i, 
region.  This  algorithm  did  give  a  better  fit  to  their  data  at  l.ig';. 
water  vapour  pressures  but  takes  no  account  of  ambient  temperature  which 
varied  over  S°C  to  30°C.  Caution  is  necessary  in  the  use  of  algorithms 
which  are  based  purely  on  field  measurements;  they  should  only  be  used 
for  conditions  under  which  the  measurements  are  made  and  not  over  the 
large  range  of  conditions  generally  experienced  in  the  atmosphere  tor 
which  they  have  not  been  tested. 

Figure  19  shows  plots  of  effective  transmittance  data  versus  visibility 
within  selected  relative  humidity  intervals  for  the  3.55  to  4.00  pm 
region.  Because  of  the  low  aerosol  extinction  coefficients  measured, 
agreement  can  be  considered  to  be  reasonable.  Recently  data  have  been 
reported ( ref .  1 5 )  in  which  LOWTKAS  3B  maritime  aerosol  extinction 
coefficients  calculated  at  3.75  um  within  the  mixed  layers  agree,  to 
within  a  factor  of  1.5  to  2,  with  the  coefficient  calculated  using 
actual  measured  aerosol  size  distributions  when  measured  visibilities 
were  of  the  order  of  10  to  23  km.  (LOWTRAS  3B  aerosol  extinction 
coefficients  were  computed  for  a  fixed  RH  of  80*,). 
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To  summarise  for  this  region,  the  model  for  HjO  continuum  presently  used 

in  LOVTRAN  5  will  need  replacing  by  a  two  component  model  once  the 
temperature  dependence  of  the  absorption  has  been  determined.  Broadband 
transmission  measurements  planned  for  the  tropical  environment  using  the 
present  equipment  will  provide  field  data  needed  for  verification  of 
updated  models  for  water  vapour  concentrations  falling  within 
15  to  25  gm/m’ . 

5.3  Variation  of  atmospheric  transmittance  in  the  1.0  to  2.5  urn  region 

In  this  spectral  region  atmospheric  transmittance  variations  are  much  more 
influenced  by  atmospheric  aerosol  absorption  variations  than  water  v.ipour 
absorption  variations.  Water  vapour  absorption  is  confined  to  spectral 

line  absorption  with  evidence  that  there  is  very  little  water  vapour 

continuum  absorption.  Analysis  of  the  data  in  this  region  has  been 

confined  to  examining  the  effects  of  aerosol  absorption  on  transmittance 
which  is  in  turn  a  function  of  aerosol  number  density,  size  distribution, 
type  of  aerosol  and  the  relative  humidity.  Data  for  this  region  were 
divided  into  onshore  and  offshore  wind  groups  and  plotted  as  a  function  of 
visibility  in  four  relative  humidity  groups.  Figures  20  to  25  reproduce 
the  data  measured  with  filters  1,2,3  and  13  over  the  two  ranges,  5,03  .and 
9.05  km,  for  onshore  winds.  The  LOWTRA.S  5  curve  in  each  plot  has  been 
determined  using  the  maritime  aerosol  model  for  a  relative  humidity 

selected  at  mid-range  in  the  interval  shown.  The  temperature  was  fixed  at 
17.5‘C  and  the  visibility  measured  at  the  time  was  used  as  an  input  to  the 
LOWTRAN  5  model. 

For  the  shorter  range  of  5.03  km  ifigures  20  to  23)  a  bias  is  evide;.', 
whereby  the  measured  values  of  transmittance  tend  to  be  below  the  priM;,.t.‘. 
values  for  all  values  of  relative  humidity.  Some  of  this  difference 
doubt  can  be  accounted  for  by  the  low  flow  rate  used  in  visibility  m.eti'rs 
allowing  heating  of  the  aerosol,  and  thereby  reducing  the  efieciiw 
humidity  and  size  of  the  aerosol  particles.  Thus,  visibilities  aciu.j'.'.y 
recorded  under  these  conditions  were  higl.er  than  were  actually  preva ;  1  i 
at  the  time. 

Finally  figure  26  provides  an  example  of  the  t ransm. i tt ance  data  me  :s'..ri-,i 
with  filter  1  over  the  5.03  km  range  for  offshore  winds.  In  this  case  th' 
rural  aerosol  model  in  LOWTRAN  5  was  used  to  compute  the  curves  and  it  is 
evident  that  discrepancies  exist.  It  is  interesting  to  note  that  more  d.it.i 
Were  recorded  for  the  higher  visibilities.  For  the  same  reasons  given 
above,  the  aerosol  extintion  coefficient  has  probably  been  underestimated 
and  this  could  account  for  some  of  the  discrepancy  seen. 

5.4  Sc int i 1 iation  effects  on  incoherent  IR  transmitted  radiation 
5.4.1  Introduction 

In  this  section  only  a  brief  descriptive  review  will  be  given  of  the 
effects  of  clear  air  turbulence  on  incoherent  propagating  beams.  The 
reader  is  referred  to  more  e.xtensive  coverages  of  the  theory  in 
references  16  and  17. 

Radiation  propagating  through  the  atmosphere  is  not  only  attenuateu  by- 
absorption  and  scattering  but  is  also  refracted  and  diffracted  by 
refractive-index  fluctuations  due  tc  small  scale  variations  la  air 
temperature.  The  resulting  fluctuations  in  intensity  of  the  transmitted 
radiation  are  known  as  scintillations.  The  size  of  the  turbulent  eddies 
present  in  the  atmosphere  is  very  important.  Those  eddies  that  are 
large  compared  with  the  receiver  diameter  tend  to  deflect  the  beam  while 
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eddies  small  with  respect  to  the  receiver  diameter  diffract  tlie 
radiation  and  cause  scintillation.  In  general,  a  combination  of  tl-.ese 
two  effects  prevails  when  long  recording  times  are  considered. 

For  the  visible  wavelengths  the  source  of  scintillation  is  known  to  be 
small  scale  fluctuations  in  the  real  component  of  the  refractive  ir.de.x 
(which  is  responsible  for  dispersion  of  the  wave  front!  caused  almost 
exclusively  by  temperature  fluctuations.  However,  scintillation  of 
visible  light  can  also  be  affected  by  strong  humidity  fluctuations  whici. 
also  affect  the  refractive  inde.x.  These  fluctuations  in  roiiactive 
index  are  expressed  in  terms  of  the  refractive  index  structure  constant 
C*..  Scintillation  is  strongest  for  near-horizontal  paths  that  are 

close  to  the  earth's  surface  (C^  ~  10  m^/*). 

The  instantaneous  random  pattern  of  sc  in  t  i  1 1  a  t  ions  ,  observed  in  a  pl.ine 
parallel  to  the  wavefront,  displays  spatial  Fourier  components  with  a 
variety  of  scales.  The  predominant  scale  size  is  the  Frcsnel-zone  size, 
/IR.  Thus  for  pathlengths  of  10  km  the  predominant  scale  size  varies 
between  10  and  32  cm  for  the  IK  wavelengths  1  .ind  10  um,  respect  ivt^  ly . 
Now  if  a  receiving  aperture  less  than  some  critical  size  (eg  --  v'TS  i 
exists  such  that  the  intensity  fluctuations  are  correlated  over  its 
area,  then  the  aperture  will  behave  as  a  point  detector.  As  the 
aperture  size  increases  the  intensity  fluctuations  become  uncorre  l.n  r,; 
over  the  aperture  and  the  intensity  is  a  spatial  average  over  some 
finite  region  of  space,  and  the  fluctuations  lend  to  decrease.  Thi.s 
effect  is  known  as  "aperture  averaging”. 

5.4.2  Scintillation  of  infrared  radiation 

The  diffraction  theory  has  been  nioder.itoly  successiul  in  predicting  tn..- 
magnitude  of  visible  scintillations  allliough  sjtur.u  ion  effects  h.ive  no; 
been  so  easily  explained.  However,  there  appe.irs  to  ho  a  scarcity  o: 
scintillation  data  for  propagating  incoherent  Ik  be.ims,  part  i  c..  i  .i  r  .  y 
involving  large  transmitted  beams  and  receiving  apertures  as  used  tne 
present  equipment.  For  IK  wave  1  i-nglhs  the  sc  int  i '.  1  it  ions  si. cal; 
decrease  for  increasing  wavelength  (in  fact  by  a  factor  ot  27  fr,',- 
1=0.63  to  10.6  pm)  if  they  arise  predominali' ly  from  refractive  in.iex 
fluctuations.  For  large  receiving  apertures,  these  scintillations  w.ll 
be  "aperture  smoothed”  thus  further  reducing  their  magnitude. 

However,  for  the  IK  region  one  has  also  to  ex.imiue  llie  effiu.ts  oi 
absorption  fluctuations  arising  from  water  vapour  pressin  .n.i 
atmospheric  pressure  fluctuations.  (The  imaginary  part  c:  tin 
refractive  index  and  the  absorption  coefficient  due  to  a  singl. 
absorption  line  are  related  by  the  expression  8  =  4nunj  where  v  is  ti.e 

wave  frequency).  A  recent  theoretical  study  ( ref .  18  )  has  sliown  that 

(i)  for  strong  humidity  fluctuations  the  refractive  index 
fluctuations  dominate  in  the  IR  and 

(li)  absorption  fluctuations  contribute  far  less  lo  the  variance  oi 
log- irtons ity  than  the  refractive  fluctuations  do. 

Thus,  for  a  point  receiver  one  would  still  expect  lo  have  a  sraal'.er 
variance  than  in  the  visible  region  except  in  water  bands  where 
absorption  is  large  and  pathlengths  are  short. 

This  now  raises  the  question  with  large  receiving  apertures;  if 
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"aperture  averaging"  reduces  the  effects  of  refraction  fluctuations  will 
absorption  fluctuations  dominate?  This,  in  effect,  is  saying  that  the 
variance  of  the  log-intensity  arising  from  absorption  fluctuations  is 
relatively  unaffected  by  aperture  size.  This  is  a  question  about  which 
little  is  known  because  of  the  limited  experimental  data  available. 

During  the  measurement  of  atmospheric  transmittance  at  the  Victor  Harbor 
site,  transmission  signals  recorded  for  wavelength  regions  between 
I  and  12  pm  (see  Table  1)  generally  exhibited  scintillation  effects  ot 
varying  magnitude.  Earlier  observations  of  scintillations  during 
transmission  measurements  over  land  were  briefly  reported  in  reference  2 
and  indicated  the  frequencies  of  the  sc int i 1 1  at  ions  fell  within 
10  to  15  Hz.  In  view  of  the  above  discussion  and  the  large  receiver 
aperture  diameter  (150  mm)  with  respect  to  the  Fre.sne  1 -zone  scale  (which 
must  imply  some  "aperture  averaging”),  scintillations  of  such  magnitudes 
appear  to  arise  from  absorption  fluctuations  which  could  be  due 
predominately  to  humidity  fluctuations.  It  should  be  noted  that  the 
scintillations  observed  on  the  recorded  transmission  signals  do  not 
totally  reflect  the  magnitude  and  frequencies  of  the  scintillations. 
This  arises  from  the  fact  that  the  phase  sensitive  detector  is  a  highly 
selective  filter  for  frequencies  away  from  the  fundamental  chopping 
frequency  (-195  Hz)  and  the  integrator  has  a  long  time  constant  (-1  s'. 
Hence,  the  recorded  scintillation  data  can  only  be  used  for  a  relative 
analysis,  that  is,  comparing  scintillation  behaviour  within  one  or  more 
measurement  periods.  As  a  general  observation  one  can  say  that  the 
magnitude  of  the  scintillation  docs  vary  substantially  from  day  to  day. 
Table  6  lists  measured  relative  fluctuations  in  the  1.48  to  2.5.  u.r. . 
3.55  to  4.0  pro  and  8.20  to  11. 7d  pm  regions  over  a  5.03  km  pat'.,  '.'cr 
scintillations  observed  'oetween  October  and  December  1979.  Typic.il.y 
the  magnitude  has  varied  by  a  factor  of  4  on  different  days.  I'ariatici.^ 
in  the  relative  peak  f  luctuat  ion.s 


Pe.ik-to-Peak  Fluctuation 


Average  Signal 


X  100 


over  a  measurement  period  on  a  particular  day  were  fairly  similar  :or 
the  throe  regions  (except  for  2/10/79)  which  appears  to  indicate  that 
the  fluctuations  are  not  sensitive  to  wavelength.  If  any  dependence 
does  exist,  the  relative  fluctuations  appear  to  be  slightly  less  for  the 
1.48  to  2.50  pm  region.  On  all  but  three  of  the  days  in  Table  b  both 
the  mean  temperature  and  the  absolute  humidity  varied  little  over  the 
measuring  period.  Thus  absorption  fluctuations  must  have  been 
responsible  for  the  greater  part  of  the  signal  fluctuations  since 
"aperture  averaging”  is  expected  to  minimise  real  refractive  index 
fluctuations . 

These  conclusions  are  substantiated  by  two  further  observations. 
Firstly  raeasurementstref .  19)  have  shown  that  the  temperature  structure 
parameter  and  hence  the  refractive  index  structure  coefficient 
varies  approximately  with  height  Z  above  the  water  surface  according  to 

Z  and  Z  ^  for  unstable  and  stable  conditions  respectively.  At 
Victor  Harbor  the  path  height  varies  over  the  water  from  about  3  to  20  m 
which  would  imply  a  fall  in  0^  by  a  factor  of  0.24  for  the  heights 
involved  in  the  unstable  situations  which  were  expected  to  prevail  at 
the  time  the  measurements  were  made.  Secondly,  signal  fluctuations  did 
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TAliLK  6.  KKLATIVK  SIGNAL  FLUCTUATIONS  NKASURED  OVER  S  km  PATIILKNGTH  FOR  3 
SPECTRAL  REGIONS  BETWEEN  OCTOBER  AND  DECEMBER  ri79 


Date 

1 .48 

to  2.50  vm 

3.55 

to  4.00  pm 

8.20 

to  11.76  ym 

Local 

T  ime 

Relative 

Fluctuation 

(%) 

Local 
T  ime 

Relative 
Fluctuat ion 
(”.) 

Loca  1 
T  ime 

j  Relative 
; Fluctuation 

1324 

3.5 

1415 

1333 

;  I 

6 .  S  1 

2/10/70 

1438 

7.0 

1530 

50.0 

1450 

b  .  5  ! 

1557 

7.0 

1654 

50.0 

IblO 

;  7.0  1 

1714 

3.5 

1721 

10.5 

1055 

2.5 

1137 

5.0 

1109 

'  6.0  1 

3/10/79 

1205 

4.0 

1245 

9.0 

1218 

'  o .  0  i 

1314 

3.0 

1319 

/ _ _2.i0 _ 1 

11/10/79 

1321 

2.5 

_ 

1411 

^  2.5 

_ 1 

1718 

_ 4 

1112 

4.5 

1148 

9.5 

1122 

^  sTo 

12/ 10/79 

1335 

11.0 

1307 

10.0 

1331 

16.0 

1305 

4.0 

1352 

5.0 

1324 

5.0 

25/10/79 

1423 

3.0 

1509 

5.0 

1439 

5.0 

1534 

4.0 

lo09 

4.7 

1  545 

5.0 

f 

1324 

13.0 

1411 

18  0 

1339 

14.0 

1450 

7 . 5 

1536 

19.0 

1504 

6 . 5 

31/10/79* 

1608 

4.0 

1712 

4.0 

lo24 

5.0 

1746 

4.0 

1837 

9.5 

1800 

3.5  j 

1904 

10.0 

2001 

9.0 

1024 

10.0 

1531 

7.5 

1619 

9.0 

1544 

10.5 

1/11/79 

163^ 

5.0 

1746 

8 . 5 

1711 

10.0 

1816 

5.0 

_ 1 

1828 

8.0 

1324 

7.0 

- 1 

1413^ 

7.5 

1338 

8.0 

1440 

7.0 

1535 

9.5 

1500 

9.0  1 

7/11/79 

1608 

7.0 

1704 

7.5 

lo23 

9.0  1 

1736 

5.5 

1825 

5.0 

1749 

7.0  ; 

1857 

4.5 

1945 

5.5 

1911 

i 

■l316~‘' 

7.5 

1403 

7.5 

r330 

7.0 

1430 

7.5 

1519 

8.0 

1448 

8.5 

12/11/79 

1551 

6.0 

1635 

8.5 

1604 

8.5 

1711 

7.5 

1759 

8.5 

1725 

!  9.5 

1831 

6.5 

1919 

8.5 

1845 

9.0 

1324 

8.0 

_ 

■QE  i  a 

8.0 

1336 

1  9.0 

1447 

5.0 

6 . 5 

1503 

;  8.5 

28/11/79 

1610 

6.5 

Hfn 

8.3 

1623 

j  10.5 

1731 

4.0 

1815 

8.0 

1743 

10.0 

1852 

4.0 

1940 

6.0 

1907 

!  9-^ 

1328 

14.0 

1415 

1342 

i  12.0 

1449 

1539 

1504 

3.0 

4/12/79* 

1610 

2.5 

1648 

1622 

,  4.0 

1740 

2.0 

1824 

1752 

!  7.0 

1856 

4.0 

1939 

1910 

1  6.0 

■ 

iHiEm 

^  6.0 

mBa 

6/12/79 

1 

6.5 

*  Refer  to  figures  27  and  28  for  actual  signal  recordings. 
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not  appear  to  change  near  sunset  when  the  atmosphere  generally  can  be 
expected  to  change  from  unstable  (free  convection  regime)  to  stable 

conditions.  In  earlier  measurements  perforaed  about  2  m  above  the 
ground,  signal  fluctuations  did  show  evidence  of  falling  as  sunset  was 
approached. 

It  is  also  noteworthy  to  recall  that  in  this  period  unusual 
scintillation  behaviour  was  observed  on  three  different  days.  On 
31/10/70  and  4/12/79  cold  fronts  passed  over  tiie  Victor  Harbor  site 

while  measurements  were  being  made  which  resulted  in  a  sliarp  drop  in 
temperature  and  an  increase  in  absolute  humidity.  Figures  27  and  26 

show  the  variation  of  temperature  and  absolute  humidity  as  the  cold 
fronts  traversed  the  site.  Actual  transmission  signals  recorded  for  six 
minutes  at  the  times  indicated  show  a  large  change  in  the  relative 
fluctuation  as  the  absolute  humidity  increases  and  temperature 
decreases,  simultaneously.  In  figure  28  the  fast  drop  in  signal 
fluctuations  for  1.48  to  2. .50  ym  region  at  1447  hours  occurred  as  the 
front  arrived  at  the  site.  During  both  these  measurement  periods,  the 
sun  was  shining. 

One  possible  explanation  for  this  phenomena  is  that  tiie  approach  of  a 
cold  front  has  reduced  the  air  temperature  to  a  value  wl.ici.  is  >,lo.-.e  i o 
the  surface  temperature  of  the  sea  water.  I:  this  were  to  occur  .t  is 
expected  that  temperature  gradients  near  the  surface  would  be  sm.i'.  1  a;.u 
combined  with  a  well  mixed  air  mass,  the  effect  of  the  i.un.ielity 

fluctuations  on  log- inlens ity  variance  could  be  substantially  reduced. 


6.  CO.VCLL'Sin.VS 

The  Victor  Harbor  measurement  site  has  provided  data  on  IK  iransmiit.une  f.,r  ■. 
coastal  site  in  a  temperate  environment  which  was  chosen  because  the  onsrio:i' 
winds  are  expected  to  provide  air  masses  which  closely  reprc.scnt  an  oper.-or.i'a:. 
maritime  environment.  These  data  have  been  classified  inlo  spectral  regions 
and  various  meteorological  groups  for  analysis.  The  main  analyses  reported 
have  been  concerned  with  the  comparisons  of  the  transmittance  data  witt; 
LOWTRAN  5  predictions  for  similar  meteorological  conditions.  It  has  neiu: 
found  that  the  data  measured  at  the  site  still  exhibit  a  degree  of  scatter 
which  is  larger  than  one  would  like  and  this,  in  part,  can  only  bo  overcome  by 
collecting  a  large  amount  of  data  for  the  wide  range  of  meteorological 
conditions  that  can  prevail. 

The  measured  transmittance  data  were  compared  with  the  LOWTRAN  5  values  in  ter, 
spectral  regions.  The  comparisons  made  in  the  three  spectral  regions  within 
the  8  to  12  vm  region  have  shown  reasonable  agreement  for  absolute  humidities 
up  to  12  gm/m’ ,  where  the  temperature  range  was  restricted  and  visibilities 
were  high.  This  conclusion  is  supported  by  the  data  measured  by  G.  .Matthews 
et  al  at  San  Nicolas  Island  off  the  West  coast  of  the  USA  (private 
communication).  For  those  spectral  regions  studied  in  the  3  to  5.5  ym  region 
the  conclusions  are  different.  The  discrepancy  which  was  reported  previously 
in  the  4.4  to  5.5  ym  region  was  found  to  remain.  Evidence  now  shows  that  this 
arises  from  the  absence  of  water  continuum  absorption  in  the  model  for  this 
region.  The  inclusion  of  a  water  continuum  absorption  component  derived  by 
AFGL  into  the  LOWTRAN  5  computer  code  has  improved  the  agreement  although  some 
discrepancy  still  remained  for  the  4.4  to  4,6  ym  region.  This  could  be  due  in 
part  to  errors  in  calibrating  the  system  for  a  narrow  spectral  region  and  ii, 
determining  very  low  transmittances  for  this  region.  The  comparisons  made  m 
the  3.5  to  4.0  ym  region  showed  more  variability,  probably  arising  from 
scatter  in  the  data  which  is  generally  more  than  in  the  other  spectral  regions 
studied.  Considering  the  accuracy  of  the  data,  the  agreement  is  reasonable. 
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On  the*  other  hand  results  for  the  1  to  2.5  ym  region  showed  a  greater 
variability.  In  this  region  the  errors  are  greater.  In  addition  to 
calibration  errors,  there  are  errors  from  (1)  measuring  visibilities  and 
(2)  using  a  rather  simple  growth  model  aerosol  in  LOWTKAN'  S  to  predict  the 
effects  of  high  humidities  on  extinction  by  the  maritime  aerosol.  So  attempt 
was  made  to  factor  out  the  aerosol  component  that  contributes  to  the  total 
effective  extinction  coefficient.  This  calculation  would  rely  on  the  accuracy 
of  the  LOWTRAN  5  model  to  predict  the  molecular  extinction  coefficient  for 
given  meteorological  and  range  conditions.  As  far  as  is  known  however,  the 
molecular  spectral  absorption  data  for  the  1  to  2.5  ym  region  is  reasonably 
accurate.  Even  so  this  calculation  must  still  introduce  some  uncertainty 
which  would  be  difficult  to  estimate  accurately.  Considerably  more  data  will 
need  to  be  obtained  in  this  region  to  check  the  influence  of  maritime  aerosol 
on  extinction  during  high  humidities. 

The  observations  reported  here  on  IR  sc  int  i  1  lat  ions  reveal  that  a  s  igr, if  leant 
level  of  scintillation  from  I  to  12  ym  is  superimposed  on  the  transmission 
signals  received  by  the  radiometer.  In  a  strictly  relative  comparison,  there 
does  not  appear  to  be  any  strong  dependence  on  wavelength,  and  the 
scint  i  1 1  at  ion®  are  not  directly  attributable  to  variations  in  the  real 
refractive  index  arising  from  temperature  fluctuations.  It  is  highly  probab'e 
that  the  scintillations  observed  arc  due  to  fluctuations  in  the  imaginary 
component  of  the  refractive  index  resulting  from  humidity  fluctuations  in  the 
atmosphere.  It  would  appear  that  if  the  air  mass  is  extremely  well  mixed  (as 
occurs  with  cold  fronts),  the  scintillations  are  considerably  reduced, 

7.  FUTURE  WORK 

A  further  measurement  programme  has  been  initiated  to  gather  more  data  on  !a 
transmittance  at  a  coastal  site,  this  time  in  the  tropical  etiviroumeni  whid. 
is  predominately  characterised  by  high  water  vapour  densities  (up  to  25  gm/in^i 
and  high  relative  humditics  (>90\).  Data  arc  particularly  required  on  ti.e 
influence  of  water  continuum  absorption  at  high  water  vapour  densities,  as 
little  experimental  data  are  available  to  compare  mode i  predictions  in  the 
3  to  5  and  8  to  12  ym  regions.  In  the  case  of  high  humidities  the  extinction 
due  to  the  maritime  aerosol  (which  is  mainly  composed  of  sea-salt) 
significantly  increased  due  to  the  growth  of  such  particles.  A  laser 
extinction  measuring  system  will  provide  absolute  transmittance  data  at  0.5j 
and  1.06  ym  wavelengths  and  provided  accurate  measurements  of  relativi 
humidity  above  90%  can  be  made,  transmittance  data  at  these  two  wavelengths 
can  be  measured  as  a  function  of  this  variable.  Such  data  will  be  valuable  in 
verifying  and  updating  models  used  to  predict  the  effects  of  relative  humidity 
on  maritime  aerosol  extinction. 
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w(X.T) 
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C  (N,T) 
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Cj,(v.T) 

‘'N 

‘'T 


Tmax 


•’ax 


NOTATION 
absorption  coefficient 

ratio  of  self  to  foreign  broadening  absorption 
coefficients 

emissivity  of  field  black*body  source 
emissivity  of  calibration  black-body  source 
wavelength  (pm) 
wavenumber  (cm  ') 

3.1416 

effective  extinction  coefficient 

slope  of  regression  line 

resultant  external  error  estimate 

imaginary  component  of  refractive  index 

coefficient  of  determination 

radiant  flu.x/unit  area  of  field  black-body 
source  at  temperature  T 

radiant  flux/unit  area  of  calibmion 
black-body  at  temperature  T 

area  of  calibration  black-body  source  aperture 

area  of  receiver  aperture 

area  of  field  black-body  source  aperture 

self  broadening  absorption  coef lie  lent 

foreign  broadening  absorption  coelficienl 

refractive  index  structure  parameter 

temperature  structure  parameter 

diameter  of  field  black-body  source  aperture 

total  lower  bound  of  relative  error  (%) 

total  upper  bound  of  relative  error  (%l 

focal  length  of  collimator  mirror 

receiver  amplifier  gain  factor 

source  radiant  intensity  (w/sr) 

constant 
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R  pathlength  (km) 

K(X)  dotector  spectral  response 

R^j^  receiver  responsivity  measured  by  linear 

regression  analysis 

R^,)^  receiver  responsivity  over  wavelength 

interval  AX(V/W) 

R„  collimator  mirror  reflectance 

M 

T  temperature  (K) 

T(X,R)  atmospheric  spectral  transmittance 

^ave  average  atmospheric  transmittance 

effective  average  atmospheric  transmittance 

Tp  receiver  filter  spectral  transmittance 

V  receiver  output  voltage 

receiver  output  voltage  for  calibration  source 
aperture  area 

integrated  radiant  flux/unit  area  from  field 
source 

integrated  radiant  flux/unit  area  from 
calibration  source 
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Figure  7.  Atmospheric  transmittance  versus  visibility  for  on-shore  winds  in 
the  1,5  to  2.5  Mm  spectral  region 
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Figure  27,  Signal  fluctuations  recorded  over  5.03  km  during  the  passage 
cold  front  across  the  measurement  site  on  51/10/79 
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Figure  28.  Signal  fluctuations  recorded  over  5.03  km  during  the  passage  of  a 
cold  front  across  the  measurement  site  on  4/12/79 
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APPENDIX  I 

EFFECTIVE  TRANSMITTANCE  MEASUREMENT  EQUATION 
I.l  Derivntion  of  fnonsur**mc‘nt  oqunLton 

This  section  briefly  describes  the  derivation  of  the  equation  which  is  used 
to  convert  the  signal  voltages  developed  by  the  f i Iter/detector 
combinations  in  the  radiometer  into  values  of  effective  transmittance 

The  irradiance,  cm  *),  at  the  radiometer  over  the  spectral  bandwidth, 

AX,  due  to  the  point  source  at  range,  R,  is 


AX 


=  '’ax 


il  .  1) 


where  is  the  source  radiant  intensity  (w  sr  *)  for  the  wavelength 

interval,  AX,  and  T  ,,(AX,R)  is  the  effective  atmospheric  transmittance 
erf 

over  the  filter  bandwidth,  AX. 


Integrating  over  the  spectral  bandpass  of  the  filter,  AX,  and  substituting 
for  gives  the  total  received  signal  as  follows 


V  = 


N  %  ■ 


if  R' 


w‘(X.  T)  Tp(X)aX 


AX 


where  w*(X,T)  =  radiant  emittance  (w  cm  ’)  from  the  field 

black-body  source  at  absolute  temperature  T, 


=  emissivity  of  the  field  black-body  sorce. 


(1.2) 


Tp(X) 

A 


K 


=  spectral  response  of  the  band-pass  filter, 

=  area  of  source  aperture, 

=  area  of  receiving  optics  aperture, 

=  obscuration  due  to  "venet ian-bl ind"  type  chopper  blades, 
and 


R 


AX 


radiometer  responsivity  (V  w  *)  averaged 
over  the  filter  bandwidth,  AX. 


w*(X,T)  for  the  field  source  is  determined  in  the  laboratory  from  surface 
temperature  'asurements  (see  reference  3).  The  radiometer  responsivity 
R^j^  is  deter  i  led  in  the  laboratory(ref .  3)  using  a  standard  black-body 

source  located  at  the  focus  of  a  collimating  mirror  (focal  length  5  ra)  and 
set  to  a  temperature  near  the  average  temperature  of  the  field  source. 
Generally  3  apertures  were  used  to  determine  the  least  squares  linear  fit 
to  the  versus  data.  Thus  using  equation  (1.2)  it  follows  that; 
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ITF' 


w(X.T)  T(13,X)  Tp(X)dX 


Vr'^  ■  ^ 


AX 


V.A. 
1  1 

1 


i=l 


where  T(13,X) 


atmospheric  transmittance  for  a  13  m  calibration 
pathlength , 

radiometer  voltage  output  for  source  aperture 

A.  . 

1 

reflectance  of  the  collimating  mirror, 
focal  length  of  the  collimating  mirror. 


Note  that  K  also  appears  in  equation  (1.3)  because  the  same  chopper  is 
used.  Combining  equations  (1.2)  and  (1.3)  results  in 


T^ff(AX.R) 


^**AX  Si  _ 1_  _1_ 

F'  •  e‘W^  •  GAs  ■  ''ax 


where 

''ax  '  /ix“(’^'T)T(13,X)Tj^.(X)dX, 

W.',  =  '  ,u‘(X,T)T  (X)dX. 

AX  AX  1 

and 


i=l 


The  gain  factor  G  for  the  amplifiers  has  been  included  in  equation  (I. A). 
In  measuring  the  choice  of  temperature  for  the  standard  black-body 

source  is  not  critical  since  for  small  temperature  changes  (eg  ±10K)  around 
980  K,  the  radiant  flux  changes  with  black-body  temperature  are  almost 
linear  for  the  spectral  region  1  to  12  um  and  the  error  arising  can  he 
neglected. 

1.2  Expression  for  calculating  effective  transmittance  from  spectral 
transmittance 

The  atmospheric  transmittance  data  are  obtained  with  a  transmissorocti  r  that 
includes  a  source,  filter  and  detector  that  have  spectral  characteristics 
over  the  wavelength  interval  AX.  Hence,  the  measured  transmittance 
T  jj.(AX,R)  is  not  the  average  of  the  spectral  transmittance  T(X,R)  for  tne 


s; 
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interval  Aa  but  is,  in  fact,  the  average  value  of  the  combination  of 
T(X,R),  the  source,  filter  and  radiometer  spectral  curves.  To  compare  the 
measured  values  of  transmittance  with  the  theoretical  spectral 
transmittance  values  derived  from,  say  LOWTRAS  5  computations,  the  source, 
filter  and  radiometer  spectral  curves  must  be  determined  and  included  in 
the  calculations.  The  equation  which  combines  all  these  curves  with  T,X.k; 
takes  the  form 


w(X.T)  TpC^)  T(X.R)  R(X3dX 

'  AX _ 

TpCX)  R(X)aX 
AX 


where  R(X)  is  the  spectral  responsivity  (Vw  '  pir.'*  )  of  the  raci ; or.e ter 
detector.  If  w(X,T),  Tp(X)  and  R(X)  are  constant  over  the  spectral  banc  iX 

(which  would  generally  be  the  case  if  AX  is  small)  then: 


W  =  df  T(X.R)aX  = 

•'AX 

Except  for  filter  7,  the  values  of  T  ,,  and  T  differ  by  no  more  than 

0,025  transmittance  for  values  between  0.2  and  0.75.  For  effective 
transmittance  calculations  made  in  this  report  using  LOWTRAS  5  data,  the 
source  spectral  response  was  assumed  to  be  a  black-body  spectral  radiance 
curve  at  T=990  K,  the  spectral  curves  of  the  individual  band-pass  filters 
were  measured  with  a  spectrometer,  and  the  detector  response  curws 
supplied  by  the  manufacturers  used  for  R(X).  The  gold  coatings  on  the 
mirrors  in  the  radiometer  and  the  "Glad-Wrap"  window  were  assumed  to  have  a 
constant  spectral  response  over  the  wavelength  region  AX. 


I 


ER1.-0265-TR 


52 


APPENDIX  II 

F.RROP  ANALYSIS  FOR  THE  MEASUREMENT  OF  ATMOSPHERIC  TRANSMITTANCE 
I  1  .  1  IntrodiicL  ion 

Tlinre  is  alw.iys  difficulty  in  determining  the  precise  mag.  itude  of  tlie 
errors  both  relative  and  absolute  associated  with  the  me.  surement  of 
atmospheric  transmittance  using  transmissometers .  The  only  thi.  g  that  c.ir. 
really  be  done  is  to  identify  the  individual  sources  of  error  (both 
operational  and  instrumental),  estimate  their  magnitudes  and  then  combine 
them  according  to  existing  laws  on  errors  which  will  give  lower  and  upper 
boundaries.  The  estimation  of  component  error  magnitudes  has  to  bo  based 
on  laboratory  measurements,  field  measurements  and  operator  experience  with 
the  instrumentation. 

I  I.  2  Error  estimates 

There  are  generally  two  methods  of  estimating  errors  associated  with 
measuring  transmittance.  They  are: 

(i)  An  external  estimate  of  the  relative  and  absolute  errors  E^  which 
are  obtained  from  a  knowledge  of  the  instrumentation  and  experiment. 

Cii)  An  internal  estimate  of  the  error  E^  based  on  .inalysis  of  t!.e 

experimental  data  itself,  eg  by  determining  the  mean  value  and  siai.d.iid 
deviation  of  a  group  of  data  from  several  transmissometers  operating  .it 
the  same  time  side  by  side. 

In  this  section,  error  analysis  will  be  devoted  to  errors  of  the  extern  il 
type  as  no  comparisons  of  the  present  equipment  have  been  made  with  other 
systems . 

11.3  External  error  estimate 

The  most  appropriate  way  to  estimate  the  external  error  for  .i 
transmissometer  system  is  to  determine,  either  empirically  by  informed 
guesses  or  by  calculation,  the  errors  associated  with  each  step  of  the 
measurement  procedure.  Those  errors  can  be  combined  according  to  the  sum 
of  the  squares  method.  If  the  individual  errors  are  independent  and 
symmetrical  then  the  resultant  error  is  given  by: 


j  =  l 


where  T  is  a  function  of  individual  terms  in  the  basic  measurement 
ef  f 

equation  each  of  which  contributes  some  error  to  the  resultant  xj 
represents  the  parameter  in  T  ,,  and  e  is  the  error  component 

Ax;  J 

(=  ) .  Thus  the  total  relative  error  (%)  in  T^^^  is  given  by 
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100 


pj  = 
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N 

^ — 
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A  t  . .  \* 
ett  ' 


‘eff  1." 
j=l 


where  T  ,,  =  f  (xi  ,Xi .  .  .x.  ) 

eff  N 

as  computed  above  gives  the  lower  bound  of  the  expected  error,  ;e  the 

relative  error  is  not  likely  to  be  less  than  the  value  grver.  b\ 
equation  (11.2).  The  maximum  error  of  the  resultant  error  will  occur  wr.er. 
all  of  the  error  components  add  up  in  the  same  direction  and  will  give  ar. 
upper  bound  to  the  overall  error  estimate,  ie 


‘MAX 


At 


eff 


Ax, 


e. 

) 


100 


j  =  l 


It  is  expected  that  this  upper  limit  on  the  total  relative  error  woulu 
occur  only  a  small  percentage  of  the  lime. 

The  basic  measurement  equation  which  provides  the  actually  meas.ir>c 

effective  transmittance  (T  ,,)  from  instrumentation  par.ameters  was  derivec 

ett 

in  Appendix  I  and  is  repeated  below. 


^^AA  Si  _rLv 
eff  -  F'e’w^^Ag  •  R^G 


Equation  (II. 4)  applies  after  correct  alignment  procedures  have  heer.  rraue 
and  the  radiometer  signal  has  been  peaked  for  a  particular  spectral  iilter. 
The  main  sources  of  error  entering  into  the  process  of  making  transmittance 
measurements  can  be  summarised  as  follows: 

(i)  Alignment  errors  during  transmission  measurements. 

(ii)  Read  out  error  induced  by  atmospheric  scintillation. 

(iii)  Gain  factor  error. 

(iv)  Calibration  errors  which  include  those  from  black-body  source 
temperature,  emissivity,  transmittance  over  short  13  m  path,  collimating 
mirror  alignment  and  other  geometric  errors. 

(v)  Calibraton  errors  in  the  field  black-body  source  which  include 
errors  from  determination  of  surface  temperature,  emissivity  and 
spectral  filter  curves. 

(vi)  Radiometer  repeatabi 1 ity  error  which  includes  errors  arising  from 
electronic  and  detector  drift  and  non-linearity. 
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To  arrive  at  the  resultant  error  equation  using  the  expression  for 

in  equation  HI. 41,  the  partial  derivatives  of  must  first  be 

(..i  M-u  1  .ited  and  (.omhuied  with  the  correspond  ing  individual  component  errors, 
if  cic/c  IS  llie  relative  error  in  the  value  of  the  eimissivity  of  ti.e 
c.a  1  1  hrat  loti  souri.e,  el;..,  iheii  i  he  resiiltaiil  error  equation  [11,21  becomes: 


Thus  the 
ideiit  i  f  ied 
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error  terms 

as  follows: 


(  I  i  .  5  1 


which  .ire 


Ac/e 

Ac'/t' 

AR  /R 
m  m 


AG/G 


relative  error  in  the  equation  omissivity  of  the 
calibration  source. 

relative  error  in  the  cmissivity  of  the  field  source. 

relative  error  of  the  reflectance  of  the  collimating 
mirror  used  during  calibration. 

rel.ltivo  error  in  the  gain  of  the  radiornei<>r  electronics. 


AD  ZD 
s  s 

=  relative 

error 

in 

the 

AF/F 

=  relative 

error 

in 

col  1 

AR/R 

=  relative 

error 

in 

the 

field  source  aperture  diameter, 
imator  mirror  focal  length, 
distance  between  receiver  and  source. 


AV/V  =  relative  error  in  the  measured  value  of  radiometer  voltage 

output  for  a  specific  f ilter/dctector  combination  during 
transmission  measurements  arising  from  atmospheric 

scintillation 

and  instrument  repeatability. 


AR*x/R*\  =  total  relative  error  in  determining  radiometer  responsivity 
AX  AX 

for  a  particular  filter/detcctor  combination  during 
cal ibrat ion . 


total  relative  error  in  the  determination  of  spectral 
radiance  from  the  calibration  black-body  source. 


AW'x/W 


AX 


total  relative  error  in  the  determination  of  spectral 
radiance  from  the  field  black-body  source. 
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Thf*  rsiimatioii  of  th**  relative  errors  AW../W..  and  AW*  /W‘,  r.i  tr.o  r.uiiant 

AX  AX  AX  Aa 

flux  frorr  the  c.u  1  ihi  at  ion  and  field  black-body  sources  i  ive  iy ,  u.is 

derive-, I  .jsMift  Lii*'  .sij’*'**'*'  root  of  tlu'  sum  of  the  s(|uari*  oi  kiu>ui.  i  uii  .  i  a  la  1 
errois  wliicli  ctiuiribuu*  to  the  resultant  error  ideiit  1 1  leii  wil!i  in»- 

evaluation  of  tlie  two  inte^jrals,  namiHy; 

ti)  .,wCX,T)TfX,13)T„(X)dX.  and 

AX  t 

(ii;  ;^j^w‘(X,T)Tj,(X)dX 

Tables  II. 1(a)  and  II. 1(b)  list  the  nature  of  the  errors,  their  i:;divi<iua; 
values  and  the  estimated  resultant  relative  error  in  calibration  and  (it-ld 
source  radiant  fluxes,  respectively.  These  two  relative  errors  show  a 

strong  wavelength  dependence  which  results  when  the  black-body  radiant  flux 
becomes  more  strongly  dependent  on  a  given  temperature  change  towards 
shorter  wavelengths. 

The  overall  magnitude  of  each  of  the  errors  which  are  independent  of 
transmittance  and  wavelength  are  listed  in  Table  11.2(a)  w'r.ile  the 
magnitude  of  those  which  ar.'  capondent  on  a  particular  f  i  Iter/uetector 
combination  are  listed  in  11. 2(b).  »fter  each  of  these  values  is 

substituted  into  equation  (11.5)  the  tot,.,  relative  error  is  calculated  ar.c 
given  in  the  form  of  square  root  of  the  suit  of  t  he  squares  in  Table  I  ’  .  i . 
In  a  similar  manner  the  maximum  error  for  the  system  can  be  c,ilc.iijt.u 
using  equation  (II. 3).  The  resultant  ernr  values  are  also  listed  in 

Table  1 1.3. 

As  most  transmittance  measurements  fell  uelween  0.1  and  0.7  the  .ibsol  .ti 
errors  for  this  range  of  have  been  calculated  for  the  three  br_au 

spectral  regions  using  the  lower  bound  relative  error  and  are  listed 
Table  II. 3. 

II. 4  Errors  in  effective  extinction  coefficient 

The  transmittance  of  the  atmosphere  over  a  particular  wavelengtli  interv.l 
AX  and  pathlength  R  can  be  expressed  in  terms  of  an  effective  value  of  the 
total  extinction  coefficient  by  fhe  following  expression: 


o 


eff 


R 


(li.n) 


where  R  is  the  pathlength  in  km. 

The  total  extinction  coefficient  is  in  fact  the  sum  of  the  aerosol 
extinction  coefficient  (which  only  varies  slowly  with  wavelength)  and  the 
effective  molecular  extinction  coefficient  over  AX  for  the  given 
pathlength. 

Now  errors  in  the  measurement  of  transmittance  will  produce  errors  in 

"eff  f ferentiating  equation  (II. 6)  with  respect  to  and  dividing  by 

■^eff 


kL-02o  S ‘TK 


5b 


The  relative  transmittance  error  results  from  calibration  errors, 
variations  in  the  source  during  measurements,  etc.  Krrors  ior  ul;;.:. 
dT  =  constant  result  from  uncertainties  in  calibration  source  aperturi*. 

gain  factor,  detector  and  electronic  drift,  etc  and  are  indopeiuicnt  o* 

T  Taking  values  of  dT  =  0.04  and  dT  ^r/T  ”  O.OhS  f'stimatod  f ro::^ 

off  eff  olfelf 

Table  1 1. 2  as  representative  of  the  lower  hounds  of  the  constant  and 
relative  transmittance  errors  for  the  present  system,  the  re  1  at  u^nsh i p 
between  errors  in  effective  extinction  coefficient  and  errors  in. 
transmittance  for  the  full  range  of  transmittance  values  is  shown  in 
figure  II.  I  .  It  is  immediately  evident  that  small  errors  in  transmi  tt 
can  lead  to  Kirge  errors  in  effective  extinction  coefficients  for  both  high 
and  low  values  of  transmittance.  As  most  of  the  transmittance  values  fail 
between  0.5  and  0.7,  the  majority  of  the  calculated  effective  extinction 
coefficient  data  would  have  proportional  errors  at  least  between  20®©  and 
40%.  Therefore  care  should  be  taken  in  interpreting  any  conclusions  from 
the  measured  extinction  coefficient  data.  The  error  analysis  in  this 
section  does  imply  that  in  making  transmittance  measurements,  pathlongth 
should  be  selected  where  transmittance  values  fall  below  60®©  if  the  error 
in  is  to  be  kept  to  a  minimum.  This  is  particularly  true  for  the  3.5 

to  4.0  ym  region  where  molecular  extinction  is  generally  small. 


Estimated  constont  error  dT-  0  0^. 


Estimated  proportional  error 


^-0  085 
T 


■Sum  of  proportlortal  and 
constant  error  terms 


Tin 

_ [_  ^ 

ln(T)  T 
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TABLE  II.l  SUMMARY  OF  ESTIMATED  RELATIVE  ERRORS  FOR  CALIBRATED  ARD  FIELD 
SOURCES 


(n)  Calibration  Source 


Type  of  Error 

Magnitude 

mm 

8  to  12  urn 

Numerical  int^'graiion 

error  in  evaluating  integral.  (A) 

0.0) 

0.01 

(’  ,  0 

Wavelength  uncertainty.  (B) 

0.01 

0.0'. 

Segl igihle 

Error  in  radiant  flux  due  to 
uncertainty  in  ave  temperature 
of  2  K.  (C) 

O.Olo 

0.01 

0.0'.)3 

Uncertainty  in  calculation  of 
T(X.13).  (D) 

0.02 

0.02 

O.Ol 

Estimated  resultant  relative 
error  =  v'(A*+B*+C’+D’) 

_ 

0.029 

_ 

0.02b 

0.011 

(b)  Field  source 


Type  of  Error 

M.ign  i  tude 

_ i 

1  to  3  ym 

3  to  5  ym 

8  to  12  ym 

Numerical  integration  error  in 
evaluating  integral.  (A) 

O.OI 

0.01 

0.01 

Wavelength  uncertainty.  (B) 

0.01 

0.01 

- 

Contribution  to  radiant 
flux  from  aperture  surround 
being  lO'C  above  ambient.  (C) 

0.005 

Error  in  radiant  flux  due  to 
uncertainty  in  ave  temperature 
of  6  K  between  calibrations.  (D) 

0.048 

0.026 

0.013 

Estimated  resultant  relative 
error  =  /(A’+B’+C’+D* ) 

0.050 

0.028 

0.017 
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TABLE  1 1. 2  SUMMARY  OF  ESTIMATED  ISDIVIDLAL  ERROR  OJMrOSESTS 
(a)  Independent  of  transmittance  and  wavelength 


(b)  Dependent  on  wavelength  and  transmittance 


Relative 

Magnitude 

1  to  3  urn 

3  to  5  um 

8  to  12  'em  I 

AV/V 

0.03 

0.03 

0.03 

0.02 

0.02 

0.02 

AR  /R 
m  m 

0.05 

0.07 

0.07 

AG/G 

0.01 

0.006 

0.U13 

*''ax''''ax’“’ 

0.029 

0.026 

0.011 

*''ax/“ax" 

0.050 

0.028 

0.017 

*Rcfer  to  Table  II. 1 
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TABLE  11.3  SUMMARY  OF  ESTIMATED  RESULTANT  TRANSMITTANCE  ERRORS 


Type  of  Error 

M.ifili  itude 

1  to  3  ym 

8  to  12  um 

IiOw»‘r  bound  rnlaiive  error 
(Square  root  of  sum  of  squares) 

0.095 

0.097 

0.093 

Ipper  bound  relative  error 

0.241 

0 .232 

0.213 

Absolute  error  (for  0.1<=  T 

eft 

<=  0.7  with  lower  bound  relative 
error) 

0.01  to  0.07 

_ 
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